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ABSTB1CT 


’"The  design,  construction,  and  testing  of  a  freestanding 
radio  frequency  (RF)  noise  measurement  and  recording  system 
in  the  30  MHz  and  150  HHz  range  is  presented.  Placement  of 
the  system  along  the  San  Andreas  fault  to  permit  the  estab¬ 
lishment  of  a  correlation  between  increased  background  RF 
noise  to  earthquake  fault  activity  for  potential  as  an 
earthquake  prediction  tool  is  described. 
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I.  INTRODUCTION 


i.  EVOLUTION  OF  EARTHQUAKE  PREDICTION 

The  problem  involving  searching  for  a  reliable  precursor 
to  earthquakes  is  being  given  much  attention  throughout  the 
world.  With  the  reality  that  within  minutes  an  earthquake 
can  trigger  an  incredibly  disasterous  act  of  nature  leaving 
death  tolls  in  the  tens  of  thousands  and  property  damage  in 
the  billions  anywhere  in  the  world.  The  desire  for 
prediction  capabilities  has  evolved  into  an  intense  investi¬ 
gative  scientific  undertaking.  Until  recently,  seismic 
activity  data  was  collected  only  to  establish  where  an 
earthquake  has  the  greatest  potential  to  occur.  Hundreds  of 
potentially  destructive  seismic  shocks  escape  unnoticed 
because  they  occur  in  remote,  sparsely  populated  areas  of 
the  earth. 

The  present  state  of  the  art  in  earthquake  prediction 
nay  be  compared  to  that  of  weather  forecasting  some  50  years 
ago.  Seismological  data  acquired  must  be  continuous  and 
detailed  from  every  area  of  the  world  that  is  subject  to 
seismic  disturbances.  Various  experimental  methods  for 
prediction  of  seismic  activity  are  widely  used  such  as  elec¬ 
tric  and  magnetic  field,  tilt,  creep,  and  strain  motion 
measurements  to  establish  a  practical  data  base  £Ref.  1]. 

Evidence  of  radio  (RF)  emission  associated  with  precur¬ 
sors  to  ea.  '  'ke  activity  is  very  sparse.  Electromagnetic 
radiation  assv  ited  with  earthquakes  appears  impulsively 
[Be£.  2],  whicl  suggests  that  these  electromagnetic  waves 
are  radiated  i  a  wide  range  of  radio  frequencies.  The  char¬ 
acteristics  oi  the  electromagnetic  radiation  depend  greatly, 
however,  on  such  factors  as  the  magnitude  of  the  earthquake. 


and  the  large  stress  converts  the  insulating  rock  to  a 
conductor  or  semiconductor.  On  reaching  the  brittle  fracture 
stress,  the  material  unloads  and  the  stress  is  reduced  and 
results  in  a  retreating  shock  wave  in  the  material.  The 
shock  wave  has  a  conductor  ahead  of  it  with  free  electrons 
that  can  fall  into  the  insulator  conduction  band  behind  the 
retreating  shock  front.  The  charges  can  accumulate  by  succe- 
sive  creep-slip  motions  to  cause  the  electronic  charge  to 
produce  an  electrical  field  beyond  the  breakdown  value  in 
the  air  or  rock  surface,  causing  an  arc  with  the  attendant 
radio  emissions. 

3.  Fracture  of  Quartz- bearing  Rocks 

The  geological  research  group  at  the  Los  Alamos 
Scientific  Laboratory  were  believed  to  be  the  first  to 
detect  radio  frequencies  when  quartz-bearing  rocks  and  other 
hard  piezoelectric  materials  were  fractured  in  laboratory 
experiments.  Using  a  standard  RF  coil  with  a  ferrite  core 
placed  about  5  to  50  CM  from  the  sample  and  a  transient 
recorder  it  was  found  that  small  cracks  in  quartz  and  tour¬ 
maline  single  crystals  produce  exponentially  decaying  RF 
transients  with  a  time  constant  of  about  10"5  seconds  as 
shown  in  Figure  2.9  .  Rock  samples  were  fractured  by 
applying  uniaxial  load  between  two  parallel  steel  anvils  or 
by  pressing  a  small  steel  ball  against  a  flat  surface  of  the 
sample.  The  signal  developed  on  the  coil  was  amplified  by  a 
wideband  amplifier  and  fed  into  one  channel  of  the  transient 
recorder.  All  measurements  were  performed  under  ambient 
conditions.  Also  the  samples  were  not  treated  in  any  way 
except  for  cutting  or  coring  [Ref.  9]. 

The  fact  that  the  fracture  related  EM  emission  was 
observed  only  in  quartz  and  tourmaline  crystals  and  in 
guartz-bearing  rocks  rules  out  triboelectr ic  effects 
(frictional  electricity)  as  the  source  of  this  emission  and 
indicates  that  it  is  related  to  piezoelectricity.  The  most 


Figure  2.8  Shock-Induced  Polarization  Signals  of  Polystyrene 
and  Plexiglas  as  a  Function  of  Shock  Pressure 

sliding  must  exist  if  the  boundary  layer  consists  of  dielec¬ 
tric  or  semiconducting  matter  and  this  electrostatic  compo¬ 
nent  is  an  appreciable  part  of  the  force  of  sliding  friction 
[fief.  8]. 

Active  fault  lines  move  by  relaxation  oscillations 
and  can  be  described  by  the  nonlinear  stick-slip  motion  that 
is  generated  when  the  static  force  of  friction  is  larger 
than  the  kinetic  frictional  value.  During  the  increasing 
stress,  the  motion  of  the  fault  proceeds  uniformly  by  creep 
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Figure  2.7  Shock-Induced  Polarization  Signals 
of  Polystyrene  Above  and  Belov  Transition  Pressure 

exhibiting  two  characteristics,  that  of  relaxation  oscilla¬ 
tion.  The  first  is  a  fixed  aaplitude  and  the  second  is  a 
period  determined  by  a  relaxation  time  for  the  building  up 
of  a  limiting  value.  The  experiments  described  the  amplitude 
of  the  jerks  was  observed  to  decrease  when  the  dielectric 
breakdown  strength  of  the  material  used  as  a  lubricant  was 
diminished,  or  when,  with  the  same  lubricant,  the  velocity 
of  sliding  was  increased.  Observations  of  contact 
electrification  show  that  an  electrostatic  force  resisting 
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Figure  2.6  Shock- Induced  Polarization  Signals 
of  Plexiglas  Above  Transition  Pressure 

signals  at  low  pressures  than  at  high  pressures  where  the 
ends  of  the  measured  signals  are  more  distorted  by  curvature 
and  obliguity  of  the  shock  front. 

2.  Contact  Electrification 

Electrification  of  insulators  and  conductors  by 
iapact  are  described  in  experiments  conducted  by  Schnurmann 
and  Warlow-Davies  of  two  metals  separated  by  an  insulating 
oil.  Slow  sliding  can  proceed  by  a  series  of  jerks 
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notion  would  result  in  a  larger  polarization  of  the  material 
by  both  mechanical  and  electrical  forces  in  aggreement  with 
the  observed  increase  in  both  Ps  and  K*  at  the  phase 


Figure  2.5  Shock- Induced  Polarization  Signals 
of  Plexiglas  Below  Transition  Pressure 


transition.  Figures  2.5,  2.6,  2.7,  and  2.8  show  that  over 
the  entire  pressure  range  of  the  measurements,  the  polariza¬ 
tion  signals  can  be  described  by  a  physical  model  involving 
three  variables,  Ps,  K*,  and 'T  .  It  should  be  noted  that 
the  conputed  curves  are  in  better  agreement  with  measured 


sliding  daring  stress  increases  and  stick-slip  motion  during 
tha  stress  release. 

1.  Insulators  and  Conductors 

The  radio  emission  theory  rests  on  the  principle 
that  an  insulator  can  be  mechanically  stressed  to  become  a 
conductor  or  semiconductor  with  subseguent  charge  migration 
from  the  conducting  rocky  matezial  to  the  insulating  rock 
under  the  equivalent  electric  field  generated  by  a  shock 
wave  during  brittle  fracture  of  the  rock.  Experimental 
evidence  from  the  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland  showed  that  shock-induced  polariza¬ 
tion  is  an  increasing  function  of  pressure  [Ref.  7]. 

At  low  pressures  shock-induced  polarization  (Ps) 
becomes  small  but  is  easily  measured  down  to  approximately 
25  Kbar.  Shock-induced  polarization  of  both  plexiglas  and 
polystyrene  was  observed  at  lover  pressures,  but  the  noise 
level  associated  with  the  experiments  prevented  accurate 
measurement.  At  high  pressures  data  indicated  transitions  Ps 
began  a  more  rapid  increase,  suggesting  transitions  to  a 
state  where  the  plexiglas  and  polystyrene  were  more  easily 
polarized  by  the  shock  wave. 

The  dielectric  constant  of  the  shocked  medium,  K* , 
was  also  an  increasing  function  of  pressure.  The  values  for 
K*  were  characterized  by  appreciable  scatter,  but  indicated 
a  gradual  increase  up  to  the  transition  pressure  where  a 
more  rapid  increase  was  noted.  The  relaxation  time,  'Y  ,  was 
also  characterized  by  appreciable  scatter,  but  remained 
between  0.5  and  1.5  usee  up  to  the  transition  pressure, 
where  a  rapid  decrease  began.  If  it  is  assumed  that  polari¬ 
zation  is  caused  by  the  motion  of  atoms  or  atomic  groups, 
then  the  sudden  decrease  in  suggest  that  motion  occurs 
more  freely  in  the  high  pressure  phase,  and  that  the 
transitions  are  to  a  more  fluid  state.  Greater  freedom  of 
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energy  radiated  over  ground  at  any  one  time  may  be  guite 
small.  This  slov  process  is  referred  to  as  a  seismic  slip  or 
creep. 

Most  earthquakes  are  generated  in  zones  where  the  huge 
plates  of  the  lithosphere,  which  make  up  the  outer  layer  of 
the  earth's  surface,  as  in  Figure  2.1  and  Figure  2.2 
[Bef.  6],  are  shearing  past  each  other.  &  conceptual  model 
of  the  faulting  process,  as  in  Figure  2.3  and  Figure  2.4, 
were  largely  developed  from  investigations  of  earthquakes 
along  the  San  Andreas  fault.  Over  the  past  137  years,  more 
than  20  severe  earthquakes  have  been  recorded  in  the  state 
of  California,  mainly  along  the  San  Andreas  fault.  The  San 
Andreas  fault  which  extends  from  the  Sulf  of  California  in 
the  South  to  Cape  Mendocino  varies  in  width  from  about  a 
meter  to  1.6  KM,  and  plunges  almost  vertically  into  the 
earth  for  15  to  20  KM.  These  are  some  of  the  reasons  why  the 
San  Andreas  qualifies  for  conducting  an  intensive,  longterm 
RF  noise  measurement  investigation  to  show  anomalous 
increases  in  the  intensity  of  electromagnetic  radiation 
precursor  to  earthquakes. 

B.  SOURCES  OF  ELECTROMAGNET IC  BADIATIOH  IH  EARTHQUAKES 

The  emissions  of  RF  are  theorized  to  be  caused  by  the 
increasing  of  mechanical  pressure  at  the  fault  causing  shear 
stresses  that  convert  the  insulating  rock  to  a  conductor  (or 
semiconductor)  by  reducing  the  energy  gap  between  the 
valence  and  conduction  bands.  During  the  stress  increase, 
microfractures  cause  shock  waves  that  accelerate  the  elec¬ 
trons  in  the  conducting  bands  to  radiate  both  acoustic  and 
electromagnetic  energy.  If  the  mechanical  forces  increase  to 
the  yield  point,  the  rock  fails  in  brittle  fracture  and 
causes  a  sudden  decrease  in  the  shear  forces.  The  mechanical 
motion  is  described  by  a  stick-slip  motion  that  has  stable 
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II.  BACRGROPND 


A.  EARTHQUAKE  GROUND  MOTION 

Equally  important  to  predicting  the  occurence  of  an 
earthquake  is  determining  which  of  the  many  vays  the  ground 
is  likely  to  shake  during  the  earthquake,  hov  strong  the 
shaking  will  be,  and  hov  long  it  will  last.  Knowledge  of 
ground  motion  can  make  it  possible  to  understand  the  charac¬ 
teristics  of  the  earthquake  source.  Most  measurements  made 
during  actual  earthquakes  at  seismological  stations  are  some 
distance  from  the  source.  The  localized  source  of  an  earth¬ 
quake  was  not  observed  for  its  significance  until  the  San 
Francisco  earthquake  of  1906. 

It  was  recognized  that  slippage  along  a  fault  in  the 
earth*s  crust  caused  the  earthquake.  The  classic  study  of 
Harry  F.  Reid  of  John's  Hopkins  University  revealed  that 
rocks  parallel  to  the  fault  had  been  strained  or  sheared. 
From  the  basis  of  his  geodetic  observations  eminated  the 
elastic  rebound  theory  of  earthquakes  which  asserts  that 
rocks  are  elastic,  and  mechanical  energy  can  be  stored  in 
them  just  as  it  is  stored  in  a  compressed  spring  [Ref.  5]. 

When  two  blocks  forming  the  opposite  sides  of  a  fault 
move  by  a  small  amount,  the  motion  elastically  strains  the 
rocks  near  the  fault.  The  frictional  bond  fails  at  its 
weakest  point  when  the  stress  becomes  large  enough.  That 
point  of  initial  rupture,  called  the  hypocenter,  may  be  near 
the  surface  or  deep  below  it.  The  seismic  waves  radiate  from 
the  hypocenter  in  all  directions  producing  the  earthquakes. 
The  point  on  the  surface  of  the  earth  above  the  hypocenter 
is  the  epicenter  of  the  earthquake.  The  rocks  rebound  over 
an  interval  of  minutes,  days  or  even  years.  The  seismic 
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density  to  produce  the  el  :tric  field  and  subsequently  the 
specific  surface  area  behaves  as  a  constant  current 
generator. 

B.  FOCUS  OF  STUDY 

The  design,  construction,  and  testing  of  a  freestanding 
radio  frequency  (SF)  noise  measurement  and  recording  system 
in  the  30  MHz  and  150  MHz  range  is  presented.  Placement  of 
the  system  along  the  San  Andreas  fault  to  permit  the  estab¬ 
lishment  of  a  correlation  between  increased  back  ground  HF 
noise  to  earthquake  fault  activity  for  potential  as  an 
earthquake  prediction  tool  is  described. 
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the  depth  of  focus,  the  geology  of  the  earthquake  region, 
the  distance  between  epicenter  and  point  of  observation,  the 
condition  of  ionospheric  propagation,  and  others.  This  is 
possibly  the  reason  why  the  intensity  of  the  electromagnetic 
radiation  is  different  for  different  earthquakes,  with 
values  varying  from  tens  or  hundreds  of  microvolts  per  meter 
up  to  tens  of  hundreds  of  volts  per  meter. 

Measurements  are  presently  being  made  to  explore  the 
possibility  of  the  existence  of  electro-magnetic  fields. 
However,  there  presently  appears  to  be  no  directed  research, 
or  a  plausible  explanation  for,  reported  radio  frequency 
emission  interference  coinciding  with  an  increase  in  earth¬ 
quake  activity.  Derr  noted  an  observation  of  noise  on  the  FM 
broadcast  band  (88-108  MHz)  cccuring  simultaneously  with 
earthquake  lights  [Ref.  3].  The  reception  of  earth-based  18 
MHz  radio  noise  by  J.W.  Warwick  of  the  High  Altitude 
Observatory,  Boulder,  Co.  at  the  time  of  the  famous  Chilean 
earthquake  still  has  been  unanswered  [Ref.  4].  Research 
personnel  involved  in  RF  interference  study,  in  interviewing 
citizen  band  and  amateur  radio  operators  in  the  Hollister, 
Ca.  area,  noted  a  consistent  report  of  "increase  in  back¬ 
ground  noise  level  preceding  earthquake  activity"  in  the 
frequency  range  of  CB  (27  MHz)  and  six  and  two  meter  amateur 
bands  (50-54  MHz  and  144-148  MHz).  Reports  typically  cited  a 
noticeable  (10  dB)  increase  in  background  noise  which 
preceded  an  earthquake  by  12  to  24  hours,  returning  to 
normal  just  before  the  tremors  were  felt  (one  hour  to  15 
minutes) .  Observations  of  the  electrostatic  field  (field 
mills)  have  reported  that  diurnal  variations  at  selected 
sites  are  accompanied  by  "high  frequency,  spiky  noise"  that 
apparently  have  no  explanation.  Various  theories  have  been 
approached  on  what  is  actually  occuring  at  the  fissure 
boundary,  whether  impact  of  dissimilar  materials  or  shock 
induced  polarization  have  developed  the  necessary  charge 
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Fig.  2.  Spectral  content  of  fracture-related 
St  emission  from  (mean  grain  size  shown  in 
brackets):  a.  coarse-grain  granite  (15  mn).;  b 
granodiorite  (2  ran);  c.  fine-grain  granite 
(from  Westerly,  R.  X.,  0.5  mm);  d.  fine-grain 
quartzite  (0.1  ram);  e.  a  2xlxC. 5  cm 
tourmaline  single  crystal  cracked  by  applying 
pressure  along  the  c  axis.  The  curves  show 
averages  of  10  spectra  for  each  sample. 


probable  mechanism  for  this  emission  is  the  rapid  drop  in 
the  piezoelectric  field  accompanying  the  sudden  release  of 
stress  when  fracture  occurs  as  shown  in  Figure  2.9  .  This 
hypothesis  is  supported  by  the  observation  that  the  ampli¬ 
tude  of  the  resulting  RF  signal  is  independent  of  load, 
since  regardless  of  the  stress  distribution  within  the 
sample  fracture  occurs  only  when  the  local  stress  exceeds  a 
specific  value  {much  higher  than  the  average  stress  in  the 
experiments)  and  therefore  drop  in  stress  is  essentially 
independent  of  load. 

Figure  2.10  compares  the  spectral  content  of  frac¬ 
ture  related  EM  emission  from  several  samples,  obtained  by 
averaging  10  different  spectra  for  each  sample.  It  is  clear 
that  a  shift  to  higher  freguency  occurs  as  the  grain  size  of 
quart  2- bearing  rocks  decrease. 

Detection  and  analysis  of  fracture-related  RF 
signals  was  very  difficult  at  distances  greater  than  50  cm. 
With  regard  to  RF  emissions  outside  of  the  laboratory,  no 
conclusions  could  be  made  from  the  study  concerning  whether 
a  substantial  increase  in  the  size  of  the  fracture  results 
also  in  a  corresponding  increase  in  the  emitted  SF  energy, 
thereby  permitting  detection  at  larger  distances. 

4 .  Porous  Rock 


Electrokinet ic  phenometa  induced  by  ground  water 
flow  associated  with  earthquakes  are  believed  to  provide  a 
possible  means  of  earthquake  prediction.  Detectable  varia¬ 
tions  of  earth  currents,  electric  potential  and  geomagnetic 
field  may  be  caused  by  diffusion  of  fluid  into  a  dilatant 
(wider  or  larger)  focal  region.  The  electrokinetic  phenomena 
controlled  by  the  processes  of  mineralized  water  migration 
through  microcapillary  rock  channels  are  an  independent 
source  of  charge  separation  and  relaxation.  The  charge 


Fig.  1.  Examples  of  EM  transients  accompanying 
fracture  of  piezoelectric  specimens  (the  digi¬ 
tal  sampling  interval  is  given  in  brackets) :  a. 
small  crack  in  a  quartz  single  crystal  (SO  ns); 
b.  small  crack  in  a  tourmaline  single  crystal 
(10  ns);  c.  failure  of  a  specimen  of  Cine  grain 
sandstone  (200  ns);  d.  failure  of  a  quartzite 
specimen,  showing  also  the  accompanying  drop  m 
load  (2030  ns,  see  text).  Note  the  different 
time  scale  for  each  transient. 


Figure  2.10  Examples  of  EH  Transients  Accompanying 
Fracture  or  Piezoelectric  Specimens 
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Figure  2. 1 1  Electric  Double  Layer  and 
Telocity  Profile  in  a  Capillary 


The  double  layer  is  Bade  up  of  a  layer  of  ions  which  are 
firmly  held  to  the  solid,  and  a  more  diffuse  mobile  layer 
extending  into  the  liquid  phase.  The  resultant  charge  of  the 
diffuse  layer  is  equal  in  magnitude  but  of  opposite  sign  to 
that  of  the  firmly  held  layer.  A  simplified  schematic 
diagram  of  the  structure  of  the  double  layer  and  profiles  of 
the  electric  potential  and  velocity  of  liquid  in  a  capillary 
is  shown  in  Figure  2.11  .  Because  of  the  electrical  charges, 
there  is  a  difference  of  electrical  potential  between  the 


solid-liquid  interface  and  the  bulk  of  the  liquid.  This  is 
called  the  electrokinetic  potential  or  zeta  potential.  When 
an  electromotive  force  (EMF)  is  applied  to  an  electric 
double  layer,  there  vill  be  a  displacement  of  the  oppositely 
charged  layers  relative  to  one  another.  If  the  solid  cannot 
move,  the  fixed  part  of  the  double  layer  also  is  unable  to 
move;  the  application  of  an  EMI  will  result  in  movement  of 
the  diffuse  layered  ions  and  liquid  [Bef.  10]. 

Applied  EMF  and  viscous  force  is  given  by: 

j  =  (£^  /q  )  grad  E  (2.1) 

where  j  is  the  fluid  flow  flux  (volume)  through  a  capillary 
per  unit  time  and  per  unit  area,  ^  is  the  zeta  potential, 
is  the  viscosity  of  the  fluid,  and  £  is  the  dielectric 
constant  of  the  fluid. 

The  streaming  potential  (E)  is  the  production  of  a 
potential  difference  when  a  liquid  is  forced  through  a 
porous  medium  or  a  capillary  tube.  It  is  given  by 

grad  E  *  -  /IjCT  )  grad  E  (2.2) 

where  and  P  are  the  electrical  conductivity  and  the  pres¬ 
sure  of  the  fluid,  respectively. 

Since  all  electrokinetic  processes  are  irreversible, 
the  thermodynamics  of  an  irreversible  process  can  be 
applied.  The  most  general  relations  between  the  electric 
current  i  and  the  fluid  flow  j,  and  the  forces  of  grad  E  and 


grad  P  are 

-i  =  111 

grad  E 

♦  LI  2 

grad  P 

(2.3) 

-j  =  L21 

grad  E 

♦  L22 

grad  P 

(2.4) 

where  Lij  are  constants.  Explicit  expressions  of  Lij’s  for  a 
capillary  tube  are: 


L 1 1  *  CT 

L12  =  121  =  -  £  ^  /  /j 

L22  =  k*  /q 


(2.5) 

(2.6) 
(2.7) 


where  k*  is  the  specific  permeability  of  the  capillary  tube. 
For  porous  media  k  is  given  by 


k  =  <^k*  (2.8) 

where  ^  is  the  porosity  of  the  porous  media.  The  basic  equa¬ 
tions  of  fluid  flow  in  the  porous  medium  are 

-I  =  cT  grad  E  -  )  grad  P  (2.9) 

-J  *  -  (461)  //J  )  grad  E  +(k/.J)  grad  P  (2.  10) 

It  is  assumed  that  the  electrical  conductivity  of  a  matrix 
substance  is  much  smaller  than  ^0~  . 

The  maximum  magnitude  of  the  current  density  associ¬ 
ated  with  the  water  flow  is  obtained  by  setting  grad  E  =  0  : 

-I  =  (4>C^  / >j  )  grad  P  (2.  11) 

This  shows  the  electric  current  is  proportional  to  the  pore 
pressure  gradient  (the  pore  pressure  is  assumed  to  change  by 
10  -  1000  bars  in  the  dilatant  focal  region  preceding  or 

following  the  earthguake)  .  If  the  mean  gradient  of  pore 
pressure  is  approximated  to  be 


grad  P  =  1  100  bar/Km 


(2. 12) 


the  electric  current  induced  is 

|I|  *  0.7  x  (10*  -  10*)  A/m*  (2.13) 

If  the  capability  of  detecting  the  water  flow  exists,  the 
direction  of  water  flow  can  be  known,  and  consequently 
predict  the  epicenter  of  the  earthquake.  The  earth  current 
due  to  electrokinetic  effects  yields  variations  of  the 
geomagnetic  field  on  the  surface. 

If  the  current  flows  uniformly  in  a  horizontal 
circular  cylinder  of  radius  R,  the  magnetic  field  caused  by 
the  current  is  simply  obtained  by  Biot  and  Savart’s  law, 

H  =  (I  /  2Tft)  •f/rR*  (2.  14) 

where  r  is  the  distance  from  the  center  of  the  cylinder  to 
the  point  at  which  the  magnetic  field  is  observed.  If 

r  *  R  *  L  /  2  (2.  15) 

then  the  field  strength  becomes 

H  *  0.2  -  200  (2.  16) 

where  the  field  strength  is  independent  of  R  and  L  for  this 
assumed  geometry.  Interestingly  enough  most  observed  magni¬ 
tudes  of  magnetic  anomalies  associated  with  earthquakes  are 
almost  all  in  the  range  estimated  above. 

C.  SOURCES  OF  ELECTROMAGNETIC  fRECURSOR  E7EHTS 

1.  Electric  Fields  and  Atmospheric  Electric  Potential 

At  the  earthquake  moment  and  just  before  it  the 
effects  due  to  disturbances  of  natural  electromagnetic  earth 


field  (NEMEF)  were  observed  repeatedly.  These  effects 
involve  quasi-stat ionary  disturbances  of  the  atmospheric 
electric  potential  (AEP)  up  to  1000  V/m,  at  distances  up  to 
some  hundred  kilometers  from  the  epicenter.  In  the  course  of 
field  observations  the  AEP  perturbations  were  noted  some 
hours  before  the  earthquake  moment  within  one  or  two  hours 
period  £Bef.  11]. 

AEP  disturbances  estimated  for  the  charged  surface 


AEP  *  2 


(2.  17) 


where 


D"  =  g  »  § 


(2.  18) 


Under  the  assumptions  that  the  linear  charge  density  q,  the 
surface  density  N,  and  the  dislocation  path  length  are: 


q  =  -  10-»°  C/m 


(2.19) 


N  =  1 0 9  -  10**  m2 


(2. 20) 


AEP  disturbances  of  100  -  1000  V/m  can  be  expected  in  the 
final  stage  of  the  earthquake  precursor  period. 

2.  Magnetic  Fields  and  Animal  Behavoir 

The  amplitude  and  time  duration  of  the  most  charac¬ 
teristic  time  varying  geomagnetic  disturbances  occurring  at 
the  earth's  surface  are  listed  in  Table  I  .  Another  way  of 
presenting  the  same  data  is  in  a  plot  of  the  spectrum  of  the 
total  magnetic  field.  Such  a  spectrum  is  shown  in  Figure 
2.12  [Bef.  12]#  where  skin  depth  {Z  =  500  (e-f)_c,'s)  is 

plotted  as  a  function  of  frequency  f  and  conductivity  xx .  It 


TABLE  I 

- 

Time  Tariations  in  the  Geomagnetic 

Field 

Characteristic 

Period 

(secs) 

Description  of 
Phenomena 

Typical 

Magnitude 

(gamma) 

10-2 

Sub-audio  frequency 
fluctuations 

0.01  sec~l 

10° 

Micro-pulsations 

0.1 

102 

Giant  pulsations 

10 

Sudden  commencements 

20 

Crochets 

20 

Bay  Phenomena 

100 

10* 

Diurnal  variation 

30 

10® 

Storm  Time  effects 

10 

Twenty-seven  day 
recurrence  effects 

8 

10® 

Annual  variation 

10 

Solar  cycle  secular 
variations 

25  year*1 

is  evident  that  the  most  doninant  signals  are  the  diurnal, 
semi-diurnal,  twenty-seven  day  (due  to  the  rotation  of  the 
sun)  and  high  frequency  changes  that  are,  in  general,  nons¬ 
tationary,  that  is  net  frequency  stable.  This  spectrum  could 
represent  also  a  crude  measure  of  discrimination  limits  for 
detection  by  animals.  It  is,  of  course,  possible  that 
animals  can  somehow  understand  these  variations  well  enough 
to  detect  changes  well  below  the  background  variations. 
Schemes  for  doing  this  are  necessary  also  for  field  experi¬ 
ments  presently  in  operation  on  the  San  Andreas  Fault  where 
information  on  the  large  scale  of  these  distrubances  is  used 
to  enhance  detection  capability  for  local  magnetic  field 
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variations  [Ref.  13].  If  animals  can  detect  effects  occur¬ 
ring  days  before  an  earthquake,  it  is  possible  that  these 
effects  involve  systematic  change  which  would  enhance  the 
spectrum  over  a  critical  frequency  range. 

There  are  a  number  of  physical  mechanisms  that  might 
be  considered  as  sources  for  earthquake-related  animal 


TABLE  II 

Order  of  Magnitude  Estimates  of  the  Various  Possible 
Physical  Hechanisms  for  Earthquake-Related  Field 

Anomalies 

A.  MAGNETO  HYDRODYNAMIC 

~  Bo 0~/uv  d  -  0.1  gamma 

(o~  =  conductivity,  yO  -  4  x  10"7,  v  =  velocity, 

D  =  length  scale/  ^  1 

B.  TELLURIC  CURRENTS 

A  ~  0" E  /  D  ™  0.01  gamma 

hr  -  conductivity,  E  =  gecelectric  field, 

D  =  depth  of  current  flow) 

C.  PIEZO  MAGNETIC  EFFECT 
AB  =  40  S  I  -  1  gamma 

(S  =  stress  change,  1=  magnetization) 

D.  LOCALIZED  CRACKS  AND  FRACTURES 
AB  =  ? 


behavior.  Table  II  lists  several  of  these  together  with 
order  of  magnitude  estimates  of  their  amplitudes.  It  would 
appear  that  for  these  effects  at  least,  it  is  difficult  to 
generate  source  amplitudes  of  more  than  a  few  tens  of 
gammas.  The  last  process  suggested  which  involves  this 
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concerns  the  relatively  unknown  electrical  and  magnetic 
effects  of  local  surface  cracks  and  fractures  might  he  a 
possible  exception  to  this. 


III.  ELECTROH AGHETIC  PROPOGATICN  CRITERIA  FOR  EARTHQUAKES 


A  geological  fault  line  becomes  a  source  of  acoustic, 
thermal  and  electromagnetic  radiations  when  the  tectonic 
plates  move  relative  to  each  ether.  The  region  near  the 
fault  line,  shown  in  Figure  3.1  is  in  both  shear  and 
compressive  stress  and  the  relative  permittivity,  <cr|  ,  is 
greater  here  than  in  the  unstressed  region  having  permit¬ 
tivity  .  The  fault  width  of  2Yo  is  the  region  of  uniform 
stress  and  the  regions  beyond  this  point  are  unstressed. 

When  fault  line  motion  commences,  a  voltage  is  generated 
causing  current  flow  across  an  interface  between  two 
different  materials.  The  interface  can  be  formed  by  the 
following  means: 

a.  Forces  near  100  Kbars  can  convert  the  insulating  rock 
to  a  semiconductor  (or  conductor)  in  solid  form  [Ref.  14]  so 
that  charges  can  flow  at  the  insulating-semiconductor  inter¬ 
face  which  is  also  the  voltage  generating  mechanism  since 
the  work  functions  differ  in  the  two  media. 

b.  Forces  above  100  Kbars  may  cause  plastic  rock  flew  to 
produce  a  liquid-solid  interface  [Ref.  15]  with  differing 
work  functions  causing  current  flow  and  charge  storage  in 
the  solid  portion  having  insulating  properties. 

c.  Shock  waves  can  distort  the  dipole  fields  in  the  rock 
to  produce  a  polarization  in  the  material  and  increase  the 
permittivity  by  nearly  an  order  of  magnitude  [Ref.  16],  The 
transient  charge  is  stored  on  the  lower  permittivity  regions 
and  accumulates  with  subsequent  shocks. 


IF  Audio  Module  Schematic 
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The  10.7  MHz  signal  is  converted  to  455  KHz  and 
amplified  by  three  stages.  Each  stage  consists  of  several 
transistors  biased  from  an  integral  zener  regulated  source. 
An  optimally  coupled  three  stage  IF  filter  provides  selec¬ 
tivity  against  adjacent  channel  interference.  Figure  4.5 
shows  a  schematic  diagram  of  the  IF  Audio  Module.  Appendix 
A  contains  Table  IV  and  Table  V,  data  of  receiver  signal 
output  voltages  for  a  given  RF  signal  input  voltage  at  both 
VHF  assigned  frequencies.  The  output  capabilities  of  the 
receiver  were  determined  from  this  information  on  the 
S-meter  (output  signal  to  the  chart  recorder). 

2 .  Recorders 

The  Rustrak  DC  recorder  prints  through  the  impinging 
action  of  its  stylus  driven  by  the  chopper  bar  against 
pressure-sensitive  chart  paper.  Its  presentation  is  a 
series  of  dots  appearing  as  a  continuous  line.  Writing 
speed  varies  with  motor  speed.  Chart  speed  and  trace  density 
depend  on  the  ratio  of  the  interchangeable  gear  box  which 
couples  the  paper  drive  to  the  motor.  Table  III  shows  the 
chart  speed  and  trace  density.  Figure  4.6  is  a  schematic  of 
the  Model  2146  dual  trace  chart  recorder  and  Figure  4.7  is  a 
schematic  of  the  same  model  with  a  signal  processor  included 
as  shown.  Both  versions  of  Model  2146  had  to  be  modified  for 
single  channel  use.  The  single  channel  modification  (S-MOD) 
involved  bypassing  a  switch  (S-2)  by  putting  a  small  jumper 
wire  between  the  normally  closed  (N/C)  and  the  normally  open 
(N/0)  terminals  as  shown  in  Figure  4.8  .  Three  of  the 
recorders  have  signal  processors  to  desensitize  the  signal 
input  to  one  volt  per  division.  The  desensitizer  board  was 
demodified  by  rearranging  the  bias  stepper  connections  at 
the  Jones  plug  making  non-bias  stepper  connections  as  shown 
in  Figure  4.9  .  There  were  no  modifications  necessary  for 
the  Model  388  chart  recorder. 


low  cost  modularized  receivers  were  procured  and  modi¬ 
fied  to  receive  at  the  assigned  frequencies.  The  output  of 
the  receiver  is  offset  to  nullify  any  ambient  noise  level 
and  is  fed  into  slow  speed  chart  recorder  to  collect  any 
detected  emissions  as  a  function  of  time.  A  block  diagram  of 
the  measurement  system  is  illustrated  in  Figure  4.2  .  The 
noise  measurement  system  can  operate  from  12  volt  DC  storage 
batteries  or  a  power  supply.  The  chart  paper  and  storage 
batteries  have  to  be  replaced  periodically.  Electromagnetic 
emission  data  collected  through  this  system  will  be  related 
to  seismic  guake  data  provided  by  the  United  States 
Geological  Survey  Office,  Palo  Alto,  Ca. 

B.  SYSTEM  COMPONEHTS 
1 .  Receivers 

The  VHF  Converter  Module  is  designed  to  amplify  and 
convert  the  specific  frequency  to  the  HF  range.  The 
converter  is  supplied  with  a  10.7  MHz  intermediate  frequency 
(IF)  transformer  for  the  mixer  output  circuit.  The  mixer 
output  circuit  consists  of  a  slug  tuned  coil,  a  capacitive 
voltage  divider,  and  an  RCA  jack  to  provide  a  50  ohm  output. 
Figures  4.3,4. 4  show  schematic  diagrams  of  the  VHF  Converter 
Modules. 

The  IF  Audio  Module  comprises  a  highly  sensitive  and 
selective  IF  amplifier,  AM  detector,  audio  amplifier,  and 
squelch  system  for  use  with  the  converter  to  make  a  VHF 
communications  receiver.  The  selectivity  of  the  unit  was 
tailored  to  operate  in  the  noise  measurement  system.  The 
unit  has  automatic  gain  control  (A GC)  of  both  IF  and  RF, 
with  the  RF  AGC  being  delayed  f cr  best  response.  An  integral 
log  detector  provides  drive  fcr  the  MS"  meter  (signal  to 
Rustrak  recorder)  output. 


17.  DESCRIPTION  OF  RESEARCH 

1.  DESIGN  GOALS 

One  of  the  most  important  problems  associated  with  using 
electromagnetic  radiation  as  a  precursor  of  earth  quakes  is 
how  to  select  the  point  of  observation.  Placement  of  the 
system  far  from  industrial  areas,  but  near  potential  earth¬ 
quake  activity  was  desired.  A  preliminary  field  survey  of 
potential  site  locations  was  conducted  under  the  direction 
of  the  United  States  Department  of  the  Interior,  Geological 
Survey.  Site  location  along  the  San  Andreas  fault  was  chosen 
because  of  it  being  known  historically  and  recently  as  an 
active  earthquake  area.  The  chosen  observation  site  is  1000 
feet  from  the  San  Andreas  Rift  Zone  (fault)  located  at  121° 
23.5'  W  and  36°  45.5*  N  which  is  three  miles  east  of  San. 
Juan  Bautista  and  seven  miles  south  of  Hollister. 

It  became  very  important  to  secure  spot  window  frequen¬ 
cies  for  the  system  so  not  to  have  any  interference  from 
transmitted  radio  waves.  The  Army  Frequency  Coordinator  for 
the  Western  United  States  identified  and  assigned  to  this 
project  VHF  bands  (38.45  MHz  and  150.75  Mhz)  .  These 
frequencies,  with  a  50  KHz  bandwidth,  have  been  allocated 
with  exclusive  use  for  a  period  of  one  year. 

For  observation  of  signal  intensity  and  direction  in  the 
VHF  range,  a  data  recording  sjstem  was  positioned  with  its 
directional  sensors  (antennas)  looking  both  up  and  down  the 
fault  as  shown  in  Figure  4.1  .  The  antennas  (eight  Yagi's) 
detect  both  horizontally  and  vertically  polarized  emissions 
at  both  reserved  frequencies.  The  omnidirectional  antennas 
(two)  are  erected  for  detection  at  both  frequencies,  and 
with  vertical  polarization.  Each  of  the  ten  antennas  are 
connected  to  a  detector  (communications  receiver) . 


Since  the  critical  angle  is  small,  the  majority  of  the 
radio  emissions  are  reflected  from  the  unstressed  region  for 
hypocenters  deep  in  the  fault.  This  wide  region  of  radio 
emission  may  account  for  the  abnormal  human  and  animal 
behavior  distant  from  the  fault  line  but  within  the  stressed 
fault  region.  When  the  field  in  the  air  space  above  the 
epicenter  exceeds  30  Kv/cm  ,  the  air  breaks  down  and 
confirms  the  reports  of  lightning  before  and  during  large 
eart hguakes. 


under  the  worst  case  of  the  hypocenter  at  the  extreme  depth, 
the  electrical  field  is  reduced  slightly  more  than  1/e,  but 
easily  measured  at  the  epicenter. 

The  values  of  permittivity  and  conductivity  vary  with 
frequency  and  the  upper  limit  cf  propogating  frequency  can 
be  estimated  from  the  graph  by  Kraus  and  Carver  £fief.  17]. 
Using  the  worst  case  of  urban  ground  at  100  megahertz,  the 
skin  depth  is  about  0.8  Km  and  if  the  hypocenter  is  halfway 
down  the  15  Km  tectonic  plates,  the  attenuation  is  e~V. 
This  limits  the  radiation  to  the  range  5  X  107  <  f  <  10®  Hz. 

C.  THE  PHOPOGATIOH  CBITEEION 

Figure  3.1  shows  the  radiation  source  at  the  hypocenter 
and  the  electric  field,  parallel  to  the  incident  plane, 
makes  an  incident  angle  ei  with  the  normal  to  the  fault.  The 
reflected  wave  is  at  an  angle  Or  with  the  normal.  The 
stressed  region  has  a  width  2Yo  and  £r(  >  •  If  the 
regions  are  insulating  and  lossless,  the  reflection  coeffi¬ 
cient  is  given  by: 


Er  /  Ei 


f  . 1  a 

Cojj  Q\  -  j  e> 


■  2  : 


Coze*,  +. 


(3.21) 


(3. 22) 


The  reflected  wave,  Er,  eguals  the  incident  wave,  Ei, 
when  the  radical  in  the  numerator  vanishes.  The  critical 
value  of  the  incident  angle  for  complete  reflection  is: 


sin  Oic  =  (^2/0.)  *5  =  (6j/6r.)  '* 


(3.23) 


This  relation  obtains  for  either  perpendicular  or 
parallel  polarization  of  the  electric  field  and  for  all  £  £ 
Oic.  Using  ^r,=  25.0  for  the  stressed  region  and  £  r,  =  5.0 
for  the  unstressed  portion,  gives 


B.  SKIM  DEPTH  AND  RADIATED  FREQUENCIES 

In  Figure  3.1,  the  fault  line  region  which  has  a 
constant  compressive  and  shear  stress  is  region  one,  and  the 
unstressed  region,  away  from  fault  line,  designated  region 
two.  The  permittivities  and  conductivities  in  the  two 


regions  are: 

(3. 13) 

D7  =n7 

(3.  14) 

(3.  15) 

II 

(3.  16) 

Where  is  the  permittivity  of  free  space  and  6^,  and£ 
are  the  relative  permittivities.  The  skin  depth,  %,  is  the 
distance  where  the  radiation  is  reduced  by  the  factor  of  1/e 
from  the  initial  value.  The  skin  depth  can  be  obtained  from 
equations: 


%  =  2<Uj6,  /C-  =  2  Co6,6/&;  (3.17) 

Using : 


£  =  5x107 

Hz,  Oj  = 

1007T  x  10*  Hz,  =10-5, 

(3.  18) 

(rt  -  25.0. 

= 

8.8  xIO'iz 

(3.19) 

II 

13.9  Km 

(3.  20) 

Since  the  San  Andreas  fault  is  about  15  Km  deep,  the 
skin  depth  is  less  than  this  value  which  indicates  that 


The  parallel  E  vector  mode,  where  the  E  field  is  in  the 
incident  plane,  can  be  obtained  from  the  vector  equations  to 
obtain  the  scalar  differential  equations: 

3Hz  /  a*  =  crEy  ♦  a  <?Ey  /  at  (3.3) 

3Ey  /  at  =  -  z  /  at  (3.4) 

Let  the  arc  discharge  have  the  spectrum  fl  <  f  <  f2  and 
=  2  f  f.  If  E  =  Eo*  ,  the  scalar  field  reduces  to: 

a£y  /  ax  -  'y2  Ey  (x,t)  =  0  (3.5) 

where 


-y2  =  j co/Jff-  - 

(3.6) 

—  XK 

Assume  a  scalar  solution  with  E  =  Eo*  e 

This  gives  the  phase  shift  and  attenuation: 

•  e 

y  =  (j  cr foe  -1) 

(3.7) 

sj  (//*>)■  5  (l-jtr-/^).* 

(3.8) 

=  j  ,5(1-j<?/  2u>e+ - ) 

(3.9) 

~  j  (//«)  /5(1- j  cr/  2^t) 

(3.  10) 

-  j  '5  +  tr-/2 

(3.  11) 

The  radiating  field  has  the  form: 


Electromagnet  Vectors  at  Fault  Hypocent 


A.  THE  E1ECTHIC  FIELD  EQUATIONS 


The  source  of  the  electromagnetic  radiations  is  the 
discharge  of  the  accumulated  electrons  that  were  formed  by 
shock  waves  or  increasing  pressures  on  the  fault.  A  contin¬ 
uous  corona  discharge  can  occur  from  the  nonuniform  fields 
that  form  on  sharp  points  of  the  rock  or  a  sudden  discharge 
can  result  when  the  field  exceeds  the  breakdown  strength  of 
the  rocks  in  uniform  or  nonunifcrm  fields.  If  the  breakdown 
occurs  in  microseconds,  radio  waves  in  the  megahertz  region 
are  radiated.  This  increasing  pressure  at  a  fault  can  be 
related  to  the  strength  and  frequency  spectrum  of  radiations 
which  can  be  correlated  to  the  acoustic  data  presently  being 
recorded  at  the  earthquake  sites. 

The  fault  hypocenter,  shown  in  Figure  3. 1  is  usually 
below  the  surface  and  is  the  scurce  of  the  electromagnetic 
radiation,  shown  as  a  point  source.  The  emissions  are  prob¬ 
ably  elliptically  polarized  but  only  the  parallel  field  is 
depicted  for  clarity.  Figure  3.2  shows  the  electric  field 
vector,  E  in  volts/meter  in  the  incident  X-Y  plane  and  the 
magnetic  field  vector,  H  in  ampere s/me ter,  normal  to  the 
incident  plane.  The  energy  propcgates  on  the  transverse  mode 
(TEM)  and  the  Poynting  vector,  P  in  watts/square  meter,  on 
the  X  axis. 

The  fault  epicenter  is  shown  in  Figure  3.1  and  is  the 
location  of  the  hypocenter  at  the  fault  surface  where  the 
electromagnetic  radiations  are  a  maximum. 

The  electromagnetic  radiations  can  be  obtained  from  the 
retarded  potential  method  but  the  objective  is  to  only  indi¬ 
cate  the  attenuation  of  the  field  with  distance  which  can  be 
found  from  the  Maxwellian  equations: 

VxH=ITE+6aE/at  (3.1) 

V  x  E  =  -/SJH  /  <jt  (3.2) 


0 


TABLE  III 

Chart  Speed  and  Trace  Density 


Model  Drive  Writing  2 

Motor  Speed  30:1 

(RPM) 


6  Gear  Train  No. 

10:1  Gear  Ratio 


2194  1  1  strike/ 

4  seconds 


3  in/hr  Chart  Speed 
76  mm/hr  Chart  Speed 
252  hrs  Paper  Duration 


1  strike/  2  in/hr 

2  seconds  50  mm/hr 

378  hrs 


Chart  Speed 
Chart  Speed 
Paper  Duration 


Trace  Density 


300  Strikes/Inch 


3 •  Interface  Amplifiers 

By  design  the  receiver  output  voltages  were  approxi¬ 
mately  between  280  mV  to  320  mV  with  RF  inputs  from  approxi¬ 
mately  10  uV  to  100  uV.  Ideally  this  40  mV  range  would  drive 
the  recorder  from  zero  to  full  scale,  but  because  of  an  op 
amp  voltage  offset  only  80  %  scale  readings  were  achieved. 
An  interface  amplifier  was  designed  to  set  full  scale 
minimum  and  maximum  needle  movement  in  the  chart  recorder. 
Calibration  of  the  meter  movement  was  accomplished  by  vari¬ 
able  resistor  settings  as  shown  in  the  schematic  of  Figure 
4.10  . 


4 .  Timer 


To  compare  seismological  disturbance  data  with  the 
noise  measurement  site  recordings  accurate  time  needed  to  be 
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Figure  4.9  Bias  Stepper  Deaodification 


indicated  on  the  chart  paper.  Because  of  a  requirement  for 
ultra-low  dissipation  and  operation  using  a  battery  source, 
a  CMOS  21  stage  counter  timer  (RCA  CD4Q45A)  was  used  in 
conjunction  with  CMOS  8  stage  presettable  synchronous  down 
counters  (RCA  CD40103B)  and  a  CMOS  dual  J-K  master-slave 
flip-flop  (RCA  CD4027B)  to  obtain  the  time.  The  2.097152  MHz 
crystal  establishes  the  one  second  clock  pulse  as  shown  in 
Figure  4.11  .  The  time  and  the  length  of  the  pulse  duration 
are  set  using  8-bit  down  counters.  When  the  timer  pulses, 
the  length  of  the  pulse  indicates  the  time.  A  fifteen  minute 
pulse  will  last  for  six  seconds.  A  one  hour  pulse  will  last 
for  twelve  seconds  and  a  twenty-four  hour  pulse  will  last 
for  eighteen  seconds.  Time  can  be  initialized  as  necessary 
by  the  reset  switch. 


V.  PERFORMANCE  OF  SYSTEM 


The  informational  parameter  of  a  possible  seismic  event 
is  variations  in  the  level  of  the  electromagnetic  radiation 
envelope  or  in  the  number  of  radiation  impulses  whose  ampli¬ 
tude  exceeds  a  pre-chosen  discrimination  threshold.  The 
noise  measurement  system  is  designed  with  a  sensitivity 
range  of  0.1  uV  to  100  uV  BF  input.  The  system  was  given  a 
"shakedown  test"  on  the  roof  of  Spanagel  Hall  at  the  Naval 
Postgraduate  School  to  see  if  it  was  sensitive  to  variations 
in  the  level  of  electromagnetic  radiation  based  on  a  24  hour 
period. 

Along  the  San  Andreas  Fault,  it  is  expected  that  varia¬ 
tions  in  signal  level  will  occur  due  to  seasonal  changes. 
For  instance,  during  autumn  and  winter  the  highest  signal 
level  falls  into  the  sunless  part  of  the  day.  Spring, 
summer,  and  the  beginning  of  autumn  have  two  maxima  during 
the  24  hour  period,  the  nightly  and  the  post-noon  ones. 
Figures  5.1,  5.2,  5.3,  5.4,  and  5.5  are  recordings  which 
show  the  noise  levels  over  a  24  hour  period. 


Figure  5.1  Signal  Level  Data  froa  Noise  Measurement  System 


Figure  5.4  Signal  Level  Data  iron  Noise  Measurement  System 
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Figure  B.5  Electromagnetic  Noise  Levels  at  Suginami 


the  observation  point  was  approximately  50  Km.  The  depth  of 
focus  was  60  Km  and  the  magnitude  was  5.0.  The  noise  lovel 
at  81  KHz  increased  to  about  12  dB  above  the  threshold  level 
about  45  minutes  before  the  main  shock.  Such  a  feature  is 
guite  similar  to  the  typical  events  presented  in  Figure  B.4. 


Figure  B.4  Electromagnetic  Noise  Levels  at  Sugadaira 


anomalous  increase  in  the  intensity.  The  amplitude  exceeded 
the  normal  level  by  more  than  15  dB.  The  normal  level  was 
approximately  8  uV.  It  can  also  be  seen  in  Figure  B.4  that 
the  intensity  dropped  sharply  back  to  the  previous  level 
exactly  at  the  moment  of  the  shock. 

Another  event  occurred  at  1247  JST  on  January  28,  1981,  as 
shown  in  Figure  B.5  £Hef.  2}.  The  location  of  the  epicenter 
was  in  the  western  Ibaragi  prefecture  on  the  northeastern 
side  of  Tokyo,  and  the  distance  between  the  epicenter  and 


Figure  B.3  Electromagnetic  Noise  Levels  at  Caucasus 
C.  JAPAN 

On  Harch  31,  1980,  an  earthquake  occurred  at  1633  Japan 
Standard  Time  (JST) .  The  magnitude  was  7,  and  the  epicenter 
was  located  at  35. 4°  n  and  135.3°  E  in  the  Kyoto  prefecture. 
The  focus  of  this  earthquake  was  very  deep,  at  approximately 
480  Km.  The  distance  between  the  location  of  the  station  and 
epicenter  was  around  250  Km. 

As  shown  in  Figure  B.4  [Ref.  2],  about  half  an  hour  before 
the  shock,  the  81  KHz  data  reception  system  recorded  an 
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APPENDIX  B 

ELECTROMAGNETIC  EMISSIONS  FROM  EARTHQUAKES 

A.  COALINGA 

An  ionospheric  perturbation  that  was  produced  by  the 
Coalinga  earthquake  of  2  May,  1983,  was  detected  by  a 
network  of  high  frequency  radio  links  in  northern 
California.  The  ionospheric  refraction  regions  of  all  five 
HF  propagation  paths,  at  distances  between  160  and  285  km 
(horizontal  range)  from  the  epicenter,  as  shown  in  Figure 

B.  1,  were  all  affected  by  a  ground-motion-induced  acoustic 
pulse,  as  shown  in  Figure  B.2,  that  propagated  to  ionos¬ 
pheric  heights.  These  observations  appear  to  be  the  first 
ionospheric  disturbances  to  be  reported  this  close  to  an 
earthquake  epicenter  [Ref.  18]. 

B.  RUSSIA 

Figure  B.3  illustrates  one  example  of  the  anomalous 
increase  in  electromagnetic  radiation  possibly  related  to 
earthquake  activity  on  September  16,  1978.  In  this  case,  the 
epicenter  of  earthquake  was  located  in  central  Iran,  and  the 
receiving  equipment  was  set  up  in  a  tunnel  50  m  below  the 
ground  surface  at  Caucasus,  USSR,  1200  Km  distant.  The 
on-set  time  of  electromagnetic  radiation  at  frequencies  of 
27  KHz  and  1.63  MHz  was  about  30  minutes  before  the  main 
shock  with  magnitude  of  7.4,  as  shown  in  Figure  B.3 
[Ref.  2]. 
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B.  RICCHflEHDATIOSS 

Data  acquired  from  the  V EF  noise  measurement  system 
should  be  taken  for  at  least  a  year  observing  background 
noise  level  variations.  Since  seismological  data  can  be 
acquired  from  the  Geological  Survey  Office,  any  significant 
earthquake  in  the  area  of  the  noise  measurement  site  should 
be  immediately  compared  with  the  recorded  data.  otherwise 

data  can  be  examined  on  a  periodic  basis. 

Although  several  studies  were  presented  concerning  the 
mechanism  by  which  rocks  emit  radiowave  energy  when  they  are 
crushed,  further  investigation  is  needed  in  this  area  and 
the  following: 

1.  The  relation  between  the  intensity  of  the  background 
noise  in  the  VHF  range  and  the  location  of  the 
epicenter.  The  antennas  of  the  noise  measurement 
system  are  looking  up  and  down  the  fault  which  makes 
location  probable. 

2.  The  relation  between  the  intensity  of  the  precursor 
VHF  radio  noise  and  the  earthquake's  magnitude  and 
depth  of  focus.  The  data  obtained  may  show  similar 
variations  in  background  noise  levels  at  38  MHz  or 
150  MHz  or  both  for  earthquakes  of  similar  magnitude 
and  depth  of  focus. 

3.  Long  term  data  collection  is  needed  to  obtain 
through  order-of-magnit ude  calculations,  the  amount 
of  VHF  radio  noise  radiated  in  these  events. 


VI.  CONCLUSIONS  AND  H ECO MM BN DATIOHS 


A.  CONCLUSIONS 

The  design,  construction,  and  testing  of  an  EF  noise 
measurement  and  recording  system  in  the  30  MHz  and  150  MHz 
range  has  been  described.  Test  results  show  that  background 
noise  levels  can  be  received  and  recorded  from  tenths  of 
microvolts  to  volts.  Theory  has  shown  that  anomalous  elec¬ 
tromagnetic  radiation  levels  from  earthquakes  can  occur. 
The  noise  measurement  system  data  will  not  only  detect  the 
presence  of  increased  background  noise  levels  at  the  desig¬ 
nated  frequencies,  but  it  may  also  give  an  indication  of  the 
direction  of  reception  relative  to  the  San  Andreas  Fault  due 
to  the  polarization  of  the  antennas. 

Besides  physically  erecting  the  antennas  along  the 
fault,  the  most  difficult  problem  was  designing  a  receiver/ 
recorder  system  that  would  be  sensitive  to  slight  changes  in 
background  noise  level  at  frequencies  near  those  of  citizen 
band  and  amateur  radio  operators.  The  signal  output  of  the 
receiver  was  sent  through  an  interface  amplifier  and  the 
chart  recorders  modified  to  accomplish  this. 

Appendix  B  presents  electromagnetic  emissions  and  other 
observations  possibly  related  to  earthquake  precusor 
activity.  Measurements  of  electromagnetic  radiation  from 
earthquakes  were  detected  by  changes  in  background  noise 
levels,  but  in  the  very  low,  low,  and  high  frequency  ranges. 
The  results  of  this  research  suggest  that  electromagnetic 
radiation  measurements  has  potential  as  an  earthquake 
prediction  tool. 
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